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ABSTRACT 

Some protoplanetary disks show evidence of inner dust cavities. Recent observations of gas and dust of these so- 
called transition disks support the hypothesis that these cavities originate from particle trapping in pressure bumps. 

We present new Atacama Large Millimeter/submillimeter Array (ALMA) continuum observations at 336 GHz of two 
transition disks, SR21 and HD 135344B. In combination with previous ALMA observations from Cycle 0 at 689 GHz, 
we compare the visibility profiles at the two frequencies and calculate the spectral index (amm). The observations of 
SR 21 show a clear shift in the visibility nulls, indicating radial variations of the inner edge of the cavity at the two 
wavelengths. Notable radial variations of the spectral index are also detected for SR 21 with values of Q:rnm~3.8 — 4.2 
in the inner region (r < 35 AU) and Q:mm~2.6 — 3.0 outside. An axisymmetric ring (which we call the ring model) or 
a ring with the addition of an azimuthal Gaussian profile, for mimicking a vortex structure (which we call the vortex 
model), is assumed for fitting the disk morphology. For SR 21, the ring model better fits the emission at 336 GHz, 
conversely the vortex model better fits the 689 GHz emission. For HD 135344B, neither a significant shift in the 
null of the visibilities nor radial variations of amm are detected. Furthermore, for HD 135344B, the vortex model fits 
both frequencies better than the ring model. However, the azimuthal extent of the vortex increases with wavelength, 
contrary to model predictions for particle trapping by anticyclonic vortices. For both disks, the azimuthal variations 
of Omm remain uncertain to confirm azimuthal trapping. The comparison of the current data with a generic model of 
dust evolution that includes planet-disk interaction suggests that particles in the outer disk of SR 21 have grown to 
millimetre sizes and have accumulated in a radial pressure bump, whereas with the current resolution there is not clear 
evidence of radial trapping in HD 135344B, although it cannot be excluded either. 

Key words, accretion, accretion disk - circumstellar matter -stars: premain-sequence-protoplanetary disk-planet for¬ 
mation 


1. Introduction 

Recent observations of transition disks (disks with inner 
dust cavities) suggest that their structures may originate 
from particle trapping in pressure bumps (e.g. van der 
Marel et al. 2013; Zhang et al. 2014). Particle accumula¬ 
tion in pressure maxima has been suggested to solve the 
problem of rapid inwards drift of particles (e.g. Weiden- 
schilling 1977; Brauer et al. 2008), implying that planetary 
embryos may form in localised environments (e.g. Klahr & 
Henning 1997; Johansen et al. 2007). Pressure bumps may 
occur because of the presence of one or multiple planets in 
the disk (e.g. Dodson-Robinson & Salyk 2011; Pinilla et al. 
2012, 2015; Zhu et al. 2012), but other phenomena, such 
as dead zones may also create pressure traps and explain 
their structures (e.g. Regaly et al. 2011; Flock et al. 2015). 
Observations of transition disks reveal that dust cavities ap¬ 
pear to be smaller at shorter wavelengths (e.g. Dong et al. 
2012; Garufi et al. 2013). This spatial segregation of small 


and mm-sized grains is a natural consequence of particle 
trapping in the planet-disk interaction scenario (e.g. Rice 
et al. 2006; de Juan Ovelar et al. 2013). 

The spectral index amm (x provides infor¬ 

mation about the particle size in protoplanetary disks (see 
Testi et al. 2014, for a review). For (sub-) micron-sized 
dust, as found in the interstellar medium, amm is expected 
to have values of > 3.5 — 4.0 (e.g. Finkbeiner et al. 1999). 
When dust grows to millimetre sizes, amm is expected to 
decrease (Draine 2006; Ricci et al. 2010). Radial increases 
in amm (on 100 AU scales) have been found for individual 
disks without cavities (e.g. Guilloteau et al. 2011; Perez et 
al. 2012), consistent with radial drift. In contrast, the in¬ 
ner region of transition disks is depleted of large grains, so 
that amm would decrease with radius from the central star 
(Pinilla et al. 2014). However, these radial variations for 
ttmm in transition disk have not been spatially resolved to 
date. 
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SR 21 



HD 135344B 



Fig. 1. ALMA observations of dust continuum emission for SR 21 (top panels) and HD 135344B (bottom panels). The left columns 
show the Band 9 (689 GHz) and Band 7 (336 GHz) continuum, with the contour overlaid at 10,20,...,90% of the peak. The third 
column shows the overlay of the contours with a step of 20% of the peak (black is Band 9, red is Band 7). The right column shows 
the resolved Omm in colour contours. The beam size is indicated in each plot. 


In this paper, we combine observations from ALMA Cy¬ 
cle 0 at 689 GHz (^450 /rm) and Cycle 1 at 336 GHz 
(^ 870^m) of SR 21 and HD 135344B (also known as 
SAO 206462) to compare the morphological structures at 
the two frequencies and calculate Omm throughout the disk. 
In addition, we compare the results with a generic model 
of particle trapping and dust growth in a pressure bump 
induced by an embedded planet, and place the results in 
the context of other transition disks. 

2. Targets and observations 

SR 21 is a G3 star located in the Ophiuchus star-forming 
region at a distance of d = 120 pc (Loinard et al. 2008). 
The disk was identified as a transition disk by its mid- 
infrared SED, with a cavity radius of ^18 AU (Brown et al. 
2007). Observations with the Submillimeter Array (SMA) 
at 850 pm confirmed a ^35 AU cavity and suggested an 
azimuthal disk asymmetry (Brown et al. 2009; Andrews et 
al. 2011), which was much more clearly seen with ALMA 
observations in Cycle 0 (Perez et al. 2014). While 450 pm 
imaging with ALMA shows strong depletion of large mm- 
sized grains, by a factor of 10^ or more in the inner region 
(van der Marel et al. 2015), H-band scattered light imag¬ 
ing indicates that the small grains are much less depleted, 
perhaps by as little as a factor of 10 (Follette et al. 2013). 

HD 135344B is an F4 star located in the Sco OB2-3 
star-forming region at a distance of d = 140 pc (van Boekel 
et al. 2005). A cavity radius of ^45 AU was also identified 
by its SED (Brown et al. 2007) and confirmed by SMA 
observations (Brown et al. 2009). One of the most intriguing 
characteristics of this disk is its spiral arms observed in 
scattered light images (Muto et al. 2012; Garufi et al. 2013). 
The observations with VLT/NACO did not show signs of a 


Table 1. Observation properties 


Target 

Band 

V 

(GHz) 

(GHz) 

-^peak 

(mJy) 

-Giotal 

(mJy) 

a 

(mJy) 

Date 

SR 21 

7 

336 

3.2 

42.4 

349 

0.08 

26/07/2014 


9 

689 

7.5 

500 

2877 

2.0 

18/07/2012 

HD 135344B 

7 

336 

3.2 

69.1 

636 

0.13 

27/07/2014 


9 

689 

7.5 

438 

3360 

2.5 

05/07/2012 


gap in small dust grains down to 28 AU radius (Garufi et 
al. 2013). 

Our observations are from ALMA Cycle 0 pro¬ 
gram 2011.0.00724.S (P.I. Perez) and Cycle 1 program 
2012.1.00158.S (P.I. van Dishoeck), taken in Band 9 and 
Band 7, respectively. The details of the calibration are de¬ 
scribed in Perez et al. (2014) and van der Marel et al. (sub¬ 
mitted). For both observation sets, the images were ob¬ 
tained using Briggs weighting with a robust parameter of 
0.5 and and the beam is set to be the same at both bands 
by the post-processing of the images, which is the best pos¬ 
sible beam size corresponding to both the Band 7 and Band 
9 observations (0.27”). The details of the images are given 
in Table 1. 

As the data are taken almost two years apart, the centre 
of the 2014 data was shifted before overlaying the images to 
compensate for the proper motion of the stars. The proper 
motion is -12,-24 mas yr“^ for SR 21 (Makarov 2007) and - 
20,-24 mas yr“^ for HD 135344B (Hpg et al. 2000). Figure 1 
shows the continuum maps and their overlay. The shape 
and contrast of the asymmetry look similar between Band 
7 and 9 for HD 135344B. In contrast, SR 21 looks more 
symmetric in Band 7. There is an azimuthal shift in the 
continuum overlay for both sources. The position of the gain 
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calibrator (the same for both observations) shifts <0.1 mas, 
meaning that the azimuthal shift is not due to calibration. 
With respect to the SMA images in 2006/2007 (Brown et 
al. 2009), the asymmetries are also shifted. 

3. Data analysis 

3.1. Visibilities and disk morphology 

The real part of the visibilities at both frequencies is shown 
in Fig. 2. These are extracted from the continuum data and 
de-projected using i = 15° and position angle P.A.=14° for 
SR 21, and i = 20° and P.A.=63°, as derived in Pontoppi- 
dan et al. (2008) and van der Marel et al. (2015). The data 
are binned by taking the mean of the available data points 
in bins of 20 kA, with a minimum of five data points per 
bin. In Fig. 2, we show the real part of the visibilities for the 
two bands. The length of the projected baseline where the 
visibilities cross zero (the null) indicates the location inner 
edge of the cavity in a given bandpass (smaller values for 
the null mean that the cavity inner edge is further out for 
the same disk properties). The nulls for SR 21 are ~220 kA 
at Band 9 and ^250 kA at Band 7. For HD 135344B the 
null is almost at the same position at both wavelengths 
(~190 kA). The imaginary part of the visibilities are pre¬ 
sented in Appendix A. 

To fit the visibility profiles and constrain the struc¬ 
ture of both disks, we used the same morphological mod¬ 
els described in Perez et al. (2014). One model is a ring¬ 
like emission described by F{r,6) = Fji 
where is the location of the peak of the ring emis¬ 
sion, Ffi the flux density at r/j, and an the ring width 
(ring model). The other model assumes a ring together 
with an azimuthal Gaussian profile to mimic a vortex 
structure (Lyra & Lin 2013). The vortex is described by 

F{r,9) = Fv with ry and 

9v being the radius and P.A. at the peak of the vortex, 
Fy the flux density at (ry, 9y), and <j 0 y,Ury is the width 
of the vortex in the radial and azimuthal direction respec¬ 
tively (vortex model). The parameters of the best-fit model 
are found by minimising using the same Markov Chain 
Monte Carlo (MCMC) approach implemented by Perez et 
al. (2014). The results are summarised in Table 2. The er¬ 
rors from the MCMC calculations are much smaller than 
the spatial uncertainty from the observations, which is typ¬ 
ically ^ 10% of the beam size (i.e. ~3 ALf for SR 21 and 
^4 for HD 135344B) for the high signal to noise of the 
data (Table 1). For the MCMC calculations, the inclina¬ 
tion and P.A. are fixed to the values found by Perez et al. 
(2014) for both disks. The model including a vortex did 
not converge to a physical solution for SR 21 in Band 7, 
because the azimuthal extension of the vortex covers more 
than 27r. In summary, the SR 21 data are best fit with a 
ring model, while the HD 135344B data are best fit with a 
vortex model (see Fig. 2). The residuals obtained by sub¬ 
tracting the best-fit model from the data in Band 7 are 
shown for both targets in Appendix B. 

In the case of SR 21, the best-fit model indicates that 
larger grains, as traced at longer wavelengths by Band 7, 
are more spatially confined than the smaller grains traced 
at shorter wavelengths by Band 9 ((Ti^~12 AU at 870 /rm 
vs. AU at 450 /rm). Furthermore, the peak of emis¬ 

sion in the ring is at larger radii in Band 7 than in Band 
9 (rH~41 AU at 870 ^m vs. AU at 450 /rm). The 



Deprojected baseline (kA) 

Fig. 2. Real part of the visibilities at 689 GHz (~450 pro., 
Band 9) and at 336 Gz (~870 ^m, Band 7) for SR 21 (top panel) 
and HD 135344B (bottom panel). The plot includes error bars 
that are of the size of the plotting points. At both wavelengths, 
the models with the best-fit parameters are over-plotted (solid 
line), which is the ring model in the case of SR 21, and the vortex 
model for HD 135344B. 


location of the inner edge of the cavity at the two frequen¬ 
cies varies from ^29 AU in Band 9 to ~35 AU in Band 7, 
which together with the wider ring, is consistent with the 
shift of the null of the visibilities (Fig. 2). For HD 135344B, 
the vortex model is consistently the best-fit model for both 
wavelengths. The inner edge of the cavity is almost constant 
in the two bands, as shown in the visibility nulls (Fig. 2). 
However, the width of the vortex in both the radial and 
azimuthal directions is significantly larger in Band 7, and 
its aspect ratio {<Je,v!'^ry) decreases from 7.1 in Band 9 to 
3.3 in Band 7. 

3.2. Spectral index 

The position-dependent spectral index (right panels of 
Fig. 1) is calculated using both bands as amm = 
ln5^/ln^^. All data with F < 5a are excluded. The 
dust emission within the cavity is detected at a high level of 
significance (S/N~320 (t and ~250cr in Band 7, and ^120a 
and ~80 ct in Band 9, for SR 21 and HD 135344B, re¬ 
spectively), so amm can be accurately computed in the 
cavity separately from the dust ring. The uncertainty of 
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Table 2. Best-fit parameters the disk morphology models of ring or ring+vortex. 


(a) Ring model 


Target 

V 


Fr 

vr 

o-R 


(GHz) 


(pm) 

(AU) 

(AU) 


336 

2.07 

0.71 

41 

12 

SR 21 

689 

1.39 

5.91 

36 

15 


336 

3.75 

0.84 

61 

19 

HD 135344B 

689 

1.18 

5.72 

60 

18 


(b) Ring+Vortex model 


Target 

V 

(GHz) 

x"" 

Fr 

(pm) 

tr 

(AU) 

(Tr 

(AU) 

Fv 

{pro) 

rv 

(AU) 

n 

0'R,V 

(AU) 

o'e,v 

(AU) 

SR 21 

689 

1.06 

4.72 

35 

14 

4.0 

46 

178 

14 

40 

HD 135344B 

336 

1.52 

0.70 

70 

14 

0.88 

43 

172 

16 

53 


689 

1.05 

5.24 

65 

16 

7.0 

42 

194 

7.0 

47 


Notes. All data of Band 9 are identical than Perez et al. (2014). The parameters for SR 21 in Band 7 of the ring+vortex model are 
omitted because of the unphysical results. The errors from the MCMC calculations are much smaller than the spatial uncertainty 
from the observations (and therefore omitted), which is ~10% of the beam size (i.e. ~3 AU for SR 21 and ~4 AU for HD 135344B). 


ttmm is calculated with error propagation from the ob¬ 
servational uncertainty, which includes the calibration un¬ 
certainties (^20% and ^10% in Band 9 and 7, respec¬ 
tively) and the noise level or rms (cr) of the observations. 
With the high signal-to-noise of the data, the uncertainty 
in the spectral index CTamm ^ given position is domi¬ 
nated by the calibration uncertainties and approximates to 
« (lni/i/i/2)~^-\/(0.22 + 0.1^) « 0.3, in both cases. How¬ 
ever, for the relative spatial changes in ammi the system¬ 
atic calibration uncertainty does not contribute since it is 
constant across the field, and changes in inside and 
outside the dust cavity have high significance with max¬ 
imum values of (Ja^^ of « 0.01. SR 21 shows significant 
radial variations of amnu with values of Q;mm~3.8 — 4.2 
within ~0.3” radius (~35 AU) and anim"^2.6 — 3.0 outside 
(Aamm ^ 120crct„„,)- This implies a lack of mm-grains in 
the inner region (r < 35 AU), and is consistent with ac¬ 
cumulation of mm-grains in a localised region in the outer 
disk. For HD 135344B, no significant radial variations for 
c^mm are found, and the values remain within a range of 
2.6 — 3.2. Although Fig. 1 shows azimuthal variations of 
Oi-m-ax for both sources, these depend considerably on how 
the two images are overlaid in the calculation of amm (see 
Appendix C), and therefore any azimuthal variation within 
a range of a„ira'^e.q2.Q — 3.2 remains uncertain. 

4. Theoretical predictions of particle trapping 

Dust trapping, in radial and azimuthal pressure bumps, de¬ 
pends on the coupling of the particles to the gas, pressure 
gradient, and disk turbulence. The dimensionless stopping 
time, Stokes number, quantifies the coupling of the particles 
and it is defined in the mid-plane (and assuming a Gaus¬ 
sian vertical profile for the gas density) as St= apsTr/2'Eg 
(with ps being the volume density of a grain, with typical 
values of ^Ig cm“^, Blum & Wurm 2008, and Eg the local 


gas surface density). Particles with St^l feel the strongest 
gas drag and therefore they move much faster to the re¬ 
gions of pressure maxima (e.g. Birnstiel et al. 2010; Pinilla 
et al. 2012). As the gas is turbulent, however, it is expected 
that the dust is turbulently mixed by the gas. The dust dif¬ 
fusion prevents the concentration of all the particles with 
St~l inside pressure traps. The drift of particles (and hence 
the trapping) is efficient for particles with St > aviso where 
ttvisc is a dimensionless number that quantifies the disk vis¬ 
cosity (Shakura & Sunyaev 1973). From the combination of 
radial drift and dust diffusion, it is expected that particles 
with St~l are more concentrated at the pressure maximum 
than St < 1. For example, the models of dust trapping by a 
vortex predict that larger grains would be more azimuthally 
concentrated in the centre of the vortex (Birnstiel et al. 
2013; Lyra & Lin 2013), as observed in IRS 48 (van der 
Marel et al. 2015) and HD 142527 (Casassus et al. 2015). 
Although turbulence can change the broadness of the dust 
concentration inside pressure bumps, it does not change the 
location of pressure maxima, i.e. the mean radial/azimuthal 
location of the dust concentration. 

Figure 3 shows the model predictions for particle trap¬ 
ping triggered by planet-disk interaction. The gas surface 
density is obtained from hydrodynamical simulations of a 
1 Mjup planet embedded in the disk at 20 AU distance from 
the star as in Pinilla et al. (2012). The initial gas surface 
density and temperature are assumed to be a power law, 
and the disk viscosity is taken to be aviso = 10“^, as found 
for typical values of disk turbulence in simulations of MRI- 
active disks (e.g. Dzyurkevich et al. 2010). The dust den¬ 
sity distributions are obtained from dust evolution models 
that include the dynamics and coagulation of dust parti¬ 
cles Birnstiel et al. (2010). When a planet opens a gap in 
the disk, a pressure bump is formed at the outer edge of 
the gap, and large particles drift to the pressure maximum 
located in this case at ~ 30 AU. The particle sizes with dif- 
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Fig. 3. Dust density distribution when 1 Mjup is embedded at 
20 AU after 1 and 5 Myr of evolution (top and bottom panel 
respectively). The lines represent particle sizes with different 
Stokes numbers (St=[10“^, 10“^, 10~^, 1]), which are propor¬ 
tional to the gas surface density profile. 

ferent Stokes number (St=[10“^, 10“^, 10~^, 1]) are shown 
to illustrate how the radial concentration becomes narrower 
for particles with higher St. The maximum grain size inside 
the trap is limited by dust fragmentation, which depends on 
Q!visc (e.g. Birnstiel et al. 2010). For higher turbulence, the 
maximum grain size that particles can reach before frag¬ 
mentation would be lower, causing the particles to easily 
escape from the bumps and then drift radially inwards. For 
more turbulent disks, a stronger pressure bump (with e.g. 
a more massive planet) would be required for trapping to 
occur (e.g. Pinilla et al. 2015). Figure 4 shows the inten¬ 
sity profile at at 450 and 870 /rm. The emission profile is 
only slightly narrower at 870 /im since the range of particles 
sizes that are traced is comparable. Both the dust density 
distributions and the intensity profiles become narrower at 
longer times of evolution. 

Because of particle trapping, the region close to the 
planet is empty of mm/cm particles, and, therefore, amm is 
expected to increase close to the location of the gap carved 
by the planet. Figure. 4 also shows the expected radial pro¬ 
files of Ofnim calculated from the models at 450 and 870 pm 
and convolved with a 2D Gaussian profile of 35 AU diam¬ 
eter. In the inner part (r < 10 AU), amm drops because of 
the presence of the mm-particles. Because there is no total 
filtration of particles at the outer edge of the gap, small 
grains still flow through the gap and grow again to mm- 
sizes in the inner region (r < 10 AU). One way to prevent 



20 40 60 80 100 


radius (AU) 

Fig. 4. Top panel: intensity radial profiles at 450 and 870 pm 
obtained from the dust density distributions after 1 and 5 Myr 
of evolution. Bottom panel: radial profile of Omm calculated from 
the intensity profiles at 450 and 870 pm, which is convolved with 
35 AU Gaussian beam. 


the presence of mm-grains is to increase the mass of the 
planet, and completely filter all particle sizes (Zhu et al. 
2012; Pinilla et al. 2012, 2015). 

At early times of dust evolution (^1 Myr), mm/cm 
grains are still distributed in the entire outer disk (r > 
30 AU) and radial variations of Omm s-re detected after con¬ 
volution. This is not the case at 5 Myr, since the concentra¬ 
tion of mm/cm particles becomes much narrower compared 
to the spatial resolution and hence the potential radial vari¬ 
ation of Onim is smeared out and is not detected. 

5. Discussion 

With the new ALMA observations of SR 21 and 
HD 135344B, it is possible to test whether radial/azimuthal 
trapping is the cause of disk structures (cavities and asym¬ 
metries). Observationally, radial trapping can be tested by 
analysing the location of the null in the real part of the 
visibilities at different wavelengths because, in the case of 
trapping, the mm-emission is expected to show a wider ring 
at shorter wavelengths. Thus, the inner edge of the ring (or 
dust cavity) should be located closer to the star at shorter 
wavelengths (Sect. 4). Additional insights can be gained by 
analysing the wavelength-dependent morphology and the 
position-dependent spectral index. 


SR 21: The current observations of SR 21 at 336 and 
689 GHz suggest that the disk morphology at longer wave- 
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lengths is better described by a ring than by a vortex. From 
observations of ^^CO, the disk mass of SR 21 has been in¬ 
ferred to be ^12 Mjup, and the average value of the gas 
surface density close to the cavity is ~80 g cm“^ (van der 
Marel et al. 2015). Analysis of recent observations of 
and C^®0 of this disk confirmed the disk mass (van der 
Marel, submitted). Assuming that at each wavelength, par¬ 
ticles with a maximum size of < 3A are traced (e.g. Draine 
2006), observations at Band 7 and 9 traced particles with 
St~2.5 X 10“^ — 5.0 X 10“^ close to the location of the 
cavity. These particles are still expected to be affected by 
radial drift and to move towards pressure maxima. 

The fact that at a longer wavelength the morphology 
is better described by a ring than a vortex is in contradic¬ 
tion with model predictions of trapping by a vortex. Thus, 
the observed azimuthal structure is likely not caused by a 
vortex. Moreover, to sustain a long-lived vortex, the radial 
width of the vortex cannot be much higher than the disk 
scale height (< 2h, Barranco & Marcus 2005) and the fact 
that the radial width of the vortex in SR 21 is much larger 
that the disk scale height (~ x 5 /i van der Marel et al. 
2015), as found from the morphology fitting of the Band 9 
data (Table 2), also disfavours the vortex scenario. 

In the radial direction, larger grains traced at 870 pm 
are more narrowly concentrated than smaller grains traced 
at 450 pm in SR 21, as demonstrated in the shift of the 
nulls of the visibilities (see Fig. 2 and Table 2). The nar¬ 
rower concentration of larger grains is in agreement with 
particle trapping by a radial pressure bump. From measur¬ 
ing the broadness of the dust concentration inside the trap 
at different wavelengths, the disk turbulence could be esti¬ 
mated (broader implying higher turbulence) if the disk tem¬ 
perature and gas surface density profile are well constraint 
from observations. Further evidence that radial trapping 
occurs in SR 21 is given by the radial changes of the spec¬ 
tral index ammi which increase outwards, as expected from 
radial trapping models (Fig. 4). The detection of radial vari¬ 
ation of ainm with the current resolution suggests that par¬ 
ticle trapping in SR 21 has occurred in < 1 Myr because 
at longer times variations of would not be detected 
(Fig. 4). The discrepancy between model predictions and 
observations for the values of amm in the inner part of the 
disk come from the fact that we do not have total filtration 
of particles in the outer edge of the gap in our model. 

Another potential explanation for variations of amm is 
high optical depth, in which case the physical temperature 
is expected to be close to the brightness temperature. To re¬ 
produce the observed differences in amm with optical depth 
effects alone, the emission must originate from a small sur¬ 
face area. Assuming that the emission is distributed in a 
homogeneous ring whose temperature is equal to the high¬ 
est brightness temperature at the peak of emission obtained 
at the two frequencies {Tg «32 K in Band 9 and ^13 K in 
Band 7), the ring must be very narrow (<3 AU in width). 
This contradicts the observed spatial extent of the emis¬ 
sion in both bands (30-40AU for SR 21 and 60-70AU for 
HD 135344B, see Table 2). Thus, the radial variations of 
Q^mm cannot be due to optically thick emission alone. Perez 
et al. (2014) demonstrated that the emission at Band 9 
is marginally optically thick at the peak of emission, so it 
can only trace some of the mass surface density, and hence 
high angular resolution observations at longer wavelength 
are necessary to confirm our predictions at the location of 
the peak. 


HD 135344B: In this case, the model with the best-fit pa¬ 
rameters for the disk morphology at the two wavelength is 
the vortex model (Table 2). The disk mass inferred from 
^^CO, ^^CO, and C^®0 is ^20 Mj^p (van der Marel et al. 
2015, and submitted), and the average value of the gas sur¬ 
face density at the location of the dust cavity is similar to 
that in SR 21 (^80 g cm“^). Therefore, the current ob¬ 
servations trace particles with similar Stokes numbers as 
SR 21. 

From the morphological models, the azimuthal width of 
the vortex increases at longer wavelengths in contradiction 
with predictions of particle trapping by vortices. Similar to 
SR 21, the radial width of the vortex is too large compared 
to the scale height of the disk. The origin of the azimuthal 
asymmetry is inconsistent with a vortex and may be re¬ 
lated to the spiral arms observed at scattered light as also 
suggested by Perez et al. (2014). In particular the bright 
spiral in the south coincides with the location of the asym¬ 
metry in the millimetre (Garufi et al. 2013; Quanz 2015, see 
also Appendix B). Indeed, the azimuthal shift of the peak 
(Fig. 1, also observed for SR 21) may be related to the phys¬ 
ical rotation of spiral arms, but higher angular resolution 
observations are needed to confirm this prediction. 

In contrast with SR 21, neither a shift of the null of 
the visibilities nor radial variations of amm are detected for 
HD 135344B, but this does not exclude radial particle trap¬ 
ping. From the model predictions (Sect. 4), the dust density 
distribution of mm-grains is more narrowly concentrated at 
the pressure maximum after 5 Myr than at 1 Myr, and ra¬ 
dial variations in amm would not be detected with the cur¬ 
rent resolution (Fig. 4). Hence, the fact that radial changes 
in amm are not detected for HD 135344B can be related to 
the fact that any instability (e.g. planet) responsible for the 
trapping formed earlier (~5 Myr ago) in HD 135344B than 
in SR 21 (^1 Myr ago). Another possibility is that trapping 
happens in more refined structures, such as spiral arms in 
gravitationally unstable disks, in which case any variations 
of amm remains unresolved (Dipierro et al. 2015). 

Besides longer evolution times, other disk and planet 
parameters, such as planet mass, turbulence, or disk tem¬ 
perature, can also affect the gap shape and thus the ra¬ 
dial concentration of mm/cm-sized particles in pressure 
bumps, leading to a different spatial distribution of small 
and large grains. High contrast asymmetries have been ob¬ 
served in other transition disks and interpreted as vortices 
(e.g. Oph IRS 48, van der Marel et al. 2013). Detection 
of vortices in transition disks may be atypical because sev¬ 
eral parameters, such as strong turbulence or feedback from 
dust to the gas, can prevent a vortex from being long-lived 
(e.g. Ataiee et al. 2013; Zhu & Stone 2014; Raettig et al. 
2015). Even in the case where a planet originally triggered 
the formation of the vortex, an eccentric orbit or the pres¬ 
ence of addition planets can also lead to its rapid destruc¬ 
tion (Ataiee & Pinilla 2015). 

The current ALMA observations of SR 21 and 
HD 135344B suggest that anti-cyclonic vortices are not the 
origin of their low contrast azimuthal asymmetries. Obser¬ 
vations at high angular resolution at longer, optically thin 
wavelengths, which provide information about the distribu¬ 
tion of larger grains, will further constrain the origin of the 
observed dust structures in transition disks. 
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Appendix A: Imaginary part of the visibilities 

Figure A.l shows the imaginary part of the visibilities at 
689 GHz (^450 ^m, Band 9) and at 336 Gz (~870 /rm, 
Band?) for SR 21 and HD 135344B. Non-zero imaginary 
visibilities indicate the presence of an azimuthal asymmetry 
(e.g. Isella et al. 2013). For HD 135344B, there are signif¬ 
icant variations from zero at both frequencies. For SR 21, 
the non-zero values are marginal in Band 7, but are signif¬ 
icant in Band 9 data. 


SR 21 



Deprojected baseline (kA) Deprojected baseline (kA) 


Fig. A.l. Imaginary part of the visibilities at 689 GHz 
(~450 nm, Band 9) and at 336 GHz (~870 /im, Band?) for 
SR 21 (top panel) and HD 135344B (bottom panel). 


Appendix B: Residual maps 

Figure B.l illustrates the residual maps for both sources 
after subtracting the best-fit models (Table 2) to the Band 7 
data (336 GHz). For SR 21 the best-fit is described by a 
ring model, whereas for HD 135344B it is a vortex model. 
The residuals for HD 135344B show a spiral structure, as 
also suggested by Perez et al. (2014) in Band 9 (689 GHz). 
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Fig. B.l. Residual maps after subtracting the best-ht models 
from Table 2, contours start at ±3 it and space by 3cr steps. 


Appendix C: Uncertainties in the azimuthal 
variations of the spectral index 

The calculation of the spectral index (amm) from the ob¬ 
servations depends considerably on how the two images are 
overlapped. To demonstrate that the apparent azimuthal 
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variations remain uncertain, we shifted the alignment of 
the Band 7 and Band 9 images of SR21 and HD135344B by 
the positional uncertainty. Figure C.l shows Omm by shift¬ 
ing 20 mas in both vertical and horizontal direction. The 
only reliable variation of amm is in the radial direction for 
SR 21, which remains significant independent of the align¬ 
ment. In HD 135344B, there is a hint of high values of Omm 
opposite to the azimuthal asymmetry, which also remains 
independent of the alignment. 
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Fig. C.l. Calculated spectral index (amm) assuming different 
alignments of the images for SR 21 (top panel) and HD 135344B 
(bottom panel). 
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